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Band offsets of semiconductor heterostructures: a hybrid density functional 
study 
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We demonstrate the accuracy of the hybrid functional HSE06 for computing band offsets of semiconductor 
alloy heterostructures. The highlight of this study is the computation of conduction band offsets with a 
reliability that has eluded standard density functional theory. A high-quality special quasirandom structure 
models an infinite random pseudobinary alloy for constructing heterostructures along the (001) growth direc- 
tion. Our excellent results for a variety of heterostructures establish HSE06's relevance to band engineering 
of high-performance electrical and optoelectronic devices. 

PACS numbers: 73.21.-b, 78.55.Cr 



Heterostructures are ubiquitous in semiconductor tech- 
nology. For instance, AlInAs/InGaAs is used for 
quantum cascade lasers and infrared photodetectors; 
InGaP/AlGaAs for high electron mobility transistors 
(HEMTs), heterojunction bipolar transistors (HBTs), 
and phototransistors; AlInAs/InP for HEMTs; In- 
GaP/GaAs for HBTs; and InGaAs/InP for single-photon 
avalanche photodiodes and HBTs. Among the most im- 
portant properties that determine the feasibility and per- 
formance of heterostructure devices are the band offsets. 
These are the discontinuities between the valence band 
maxima (VBM) or conduction bands minima (CBM) of 
each semiconductor at their common interface, and act as 
barriers to electrical transport across the interface. Band 
engineering of novel devices with desired properties, par- 
ticularly quantum cascade lasers and quantum dot-based 
devices, critically require a precise knowledge of band off- 
sets. However, reliable measurements and predictions of 
band offsets continue to be challenging despite extensive 
theoretical and experimental efforts^— 

Density functional theory (DFT) is an efficient method 
for calculating electronic structure. The accuracy of 
DFT calculations is controlled by the choice of exchange- 
correlation (XC) functional. Local and semi- local func- 
tionals such as LDA (local density approximation) and 
PBE (Pcrdcw-Burkc-Ernzcrhof) 4 fail to produce accu- 
rate bandgaps, and in extreme cases predict small gap 
semiconductors as metals. Hybrid XC functionals, that 
include a fraction of Hartree-Fock (HF) exchange, pro- 
vide a promising alternative. In this letter, we deter- 
mine the suitability of a hybrid functional HSE06 (Heyd- 
Scuseria-Ernzerhof )& to compute band offsets of several 
III-V compounds and pseudobinary alloy heterostruc- 
tures. HSE06 includes a fraction, a, of screened, short- 
range HF exchange to improve the derivative discontinu- 
ity of the Kohn-Sham potential for integer electron num- 
bers (default HSE06 uses a=0.25). This functional was 
recently used to predict the band alignments throughout 
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FIG. 1. Band alignments for Alo.5Gao.5As/GaAs heterostruc- 
ture computed with mod-HSE06 (a = 0.30, see text) and 
PBE, in comparison with experiment.— The direct and indi- 
rect bandgaps are shown for GaAs and Alo.5Gao.5As, respec- 
tively. The hybrid functional shows significant improvement 
over PBE for bulk bandgaps and both valence and conduction 
band offsets of the heterostructure. 



the composition range of InGaN. 6 However, the bands 
were constructed and aligned with respect to vacuum in- 
stead of the conventional method of computing the band 
alignment involving a heterostructure supercell. 

Figure [TJ highlights the success of HSE06 in comput- 
ing valence and conduction band offsets of the classic 
Alo.5Gao.5As/GaAs heterostructure in close agreement 
with experiment. HSE06 also shows significant improve- 
ment over PBE for computing accurate lattice constants, 
bandgaps, and cation outermost d-orbital binding ener- 
gies for the bulk III-V compound semiconductors. 7 These 
results signal advantages of HSE06 over traditional func- 
tionals for computing electronic properties. 

Since the percentage of HF exchange in a hybrid func- 
tional is not a universal constant and the optimal value 
may be system-dependent, it is worthwhile to study the 
variation in bandgaps as a function of a in HSE06. Figure 
2 demonstrates close agreement between the computed 
and experimental direct bandgaps for III-V phosphides 
and arsenides using the default HSE06 functional. Since 
a = 0.30 describes the bandgaps of both AlAs and GaAs 
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FIG. 2. (Color online) Calculated direct bandgaps (open 
symbols) vs. fraction of Hartree-Fock mixing (a) in HSE06 
functional for III-V (a) arsenides and (b) phosphides. The 
vertical line passes through a=0.25, the default fraction for 
HSE06. PBE results are shown at a—0. Filled (partially 
filled) red symbols indicate direct (indirect) experimental 
bandgaps. The experimental bandgaps do not vary with a 
but are positioned according to the a needed in HSE06 to re- 
produce those values. Default HSE06 gives bandgaps close to 
experimental values for InAs, GaP and InP. An optimal value 
of a=0.3 is required to obtain the experimental gaps for both 
AlAs and GaAs. We use default HSE06 for A1P as different 
a values of 0.33 and 0.12 are required to match experimental 
indirect and direct bandgaps, respectively. 



TABLE I. Comparison of HSE06 and PBE calculated valence 
band offsets AE V (eV) and conduction band offsets A_E C (eV) 
of III-V binary and pseudobinary alloy heterostructures with 
experiment values. All the alloys in the present study are of 
the form A0.5B0.5C. The band offsets computed with HSE06 
show much better agreement with experimental values than 
PBE. Asterisks indicate that mod-HSE06 (a=0.30) is used 
instead of the default HSE06 functional (a=0.25). A posi- 
tive value of AE V (AE C ) for a heterostructure S1/S2 implies 
that the valence (conduction) band edge of semiconductor S2 
lies higher than that of the semiconductor SI. The valence 
and conduction band offsets have opposite signs for type-I 
(straddling) heterostructures and the same sign for type-II 
(staggered) heterostructures. 



Heterostructure 


AE V (eV) 


AB C (eV) 


S1/S2 


HSE06 PBE Exp. 


HSE06 PBE Exp. 


AlAs/GaAs 


0.52* 0.45 0.53 a 


-1.02* -1.08 -1.05 a 


AlP/GaP^ 


0.54 0.49 0.55 c 


0.58 0.57 0.38 c 


AlSb/GaSb> 


0.38 0.35 0.38^ 


-0.67 -1.21 -0.51^ 



AlGaAs/GaAs^ 
InGaP/GaAs 
InP/InGaAs 
InP/AlInAs 



0.26* 0.21 0.27^ 

0.32 0.26 0.31^ 

0.36 0.27 0.34^ 

0.16 0.14 0.17 f 



-0.42* -1.12 -0.31* 

-0.24 -0.34 -0.18^ 

-0.38 -0.42 -0.27^ 

0.22 0.17 0.25 f 



AlInP/InGaP^ 
AlInAs/InGaAs 
InGaP/AlGaAs 



0.22 0.19 0.24£ 
0.23 0.18 0.22^ 
0.11 0.08 0.09' 



-0.23 -0.74 -0.26^, 
-0.57 -0.54 -0.5V1 
0.25 0.19 0.28 j 



d/o InGaP^ 



0.01 0.02 



-0.24 -0.18 0.15 



a From Ref, 2 and references therein 

b The indirect bandgaps of A1P, GaP, AlSb, AlGaAs, AllnP arc 

used to get AE C for these heterostructures. 
c From Ref JL 

d No exp. data; AE C calculated from exp. AE V and AE g 
c From Refl. 
f From Ref. 10 . 
s From Refii. 

h From Ref 3 and references therein 

1 No exp. data; AE V calculated from exp. AE C and AE g 
J From Refi 2 .. 

d=disordered and o=ordered with Cu-Pt (Lli) ordering 
1 From Refi 3 .. 



band offset, AE V , is calculated as 



so well, a reasonable choice for a in the AlGaAs alloy 
would be 0.30, which we refer to as mod-HSE06. For all 
other systems we use the default a — 0.25 as no other 
value will work for both components of the alloy. 

We employ the average electrostatic potential tech- 
nique^ to compute the band offsets of the heterostruc- 
tures S1/S2, where SI and S2 are the semiconductors 
constructing a heterostructure. The bulk valence band 
edges are aligned through a reference potential calculated 
across the interface of the heterostructure. The difference 
in valence band maxima of S2 and SI is AiSyeM- The 
discontinuity in this reference potential across the het- 
erostructure interface is defined as AT4top- The valence 
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The conduction band offset is determined from AE V and 
the difference in bulk bandgaps, AE g , as 



A £(S2-si) = A£( S2 - S1 ) + A£ (S2 - S1) 



(2) 



TablelH by comparing the DFT and experimental band 
offsets for several III-V compounds and pseudobinary 
alloy heterostructures, establishes that both conduction 
and valence band offsets calculated by HSE06 are much 
closer to experiments than PBE. HSE06 predicts the ac- 
curate nature and magnitude of the bandgaps^ and hence 
band offsets, whereas for PBE, the error in bandgaps 



translates to error in conduction band offsets. Sam- 
pling Table H we start with AlSb/GaSb. PBE underesti- 
mates the bandgaps of both constituents,- 7 predicts GaSb 
metallic, and produces a giant AE C . HSE06 corrects the 
bandgaps as well as offsets. A major success of HSE06 
is evident for the more complex alloy heterostructures 
such as AlInP/InGaP; HSE06 predicts AllnP as indirect 
bandgap material, in agreement with experiment, and 
computes accurate offsets. PBE, on the other hand, pre- 
dicts AllnP as a direct gap semiconductor with a bandgap 
of 1.52 eV (and indirect gap of 1.72 eV). 7 The PBE cal- 
culated AE C , using the direct (indirect) bandgap, is -0.54 
eV (-0.74 eV), much bigger than the experimental value. 
Another specific example is the semiconductor InGaAs, 
which PBE predicts as metallic, 7 Although both PBE 
and HSE06 give numerically similar band offsets for A1I- 
nAs/InGaAs and InP/InGaAs, PBE fails on the physics. 
A recent study on pseudobinary alloys shows that HSE06 
accurately predicts the alloy concentration for the direct- 
indirect bandgap crossovers. 15 These results nail HSE06's 
predictive power for optoelectronics. 

We use the planewave projector augmented-wave 
(PAW) method^ with PBE and the HSE06 hybrid func- 
tional in the VASP codei^— The outermost d electrons of 
cations are treated as valence. We use a planewave energy 
cut-off of 500 eV for all our calculations. The Brillouin 
zone integration for the bulk III-V binaries and their het- 
erostructures is performed on a T-centered 8x8x8 and 
8x8x1 k-point meshes, respectively. For the pseudobi- 
nary alloys and alloy heterostructures, 4x8x4 and 4x8x2 
T-centered k-point meshes are used, respectively. Atomic 
relaxations are not taken into account for the common- 
anion materials. However, when the anions are different 
in the two semiconductors making a heterostructure, we 
relax the interfacial atoms. 

The computation of an infinite random alloy requires 
a large supercell and is arduous, particularly with the 
hybrid functional. In order to simulate such an alloy us- 
ing a finite supercell with reasonable computational ef- 
fort, we employ a special quasirandom structure (SQS), 20 
generated using the Alloy Theoretic Automated Toolkit 
(ATAT). 21 Only the 16 cations on a fee sublattice were 
distributed according to the SQS construction, whereas 
the 16 anions are located on the separate sublattice that 
makes up the 32 atom zincblende supercelli 7 - The degree 
to which this SQS matches an infinite perfect random al- 
loy is based on the behavior of the first few radial correla- 
tion functions. We search all possible 16-atom fee super- 
cells with two lattice- vectors orthogonal to (001) for easy 
construction of a heterostructure in a (001) growth direc- 
tion yielding the smallest deviation in the radial correla- 
tion functions from the perfectly random alloy. The SQS 
employed in this work has radial correlation functions 
that match the perfect random alloy up to the fourth 
nearest neighbor pairs. The substrate lattice constant, 
computed with either functional, is used for alloy het- 
erostructures. The heterostructures of alloys (binary III- 
V's) are modeled as 4+4-layer thick supercells, 64 atoms 



(16 atoms) with a (001) interface. For either functional, 
the average of the lattice constants of the nearly lattice- 
matched bulk materials, computed with that particular 
functional, is used for the heterostructure. 

To conclude, we have accurately computed the valence 
and essentially the conduction band offsets, using the hy- 
brid functional HSE06, for several heterostructures of di- 
rect technological significance. The importance of this 
study lies in its overcoming the limitations of theory 
for excited states that restricted its application to semi- 
conductor electronics. Our excellent results for a broad 
range of materials indicate that theory can predict struc- 
tures with desirable electrical and optoelectronic proper- 
ties. More resources may conquer strained and lattice- 
mismatched interfaces increasingly being used for ad- 
vanced devices. 
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